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We study the spin-down of a neutron star during its early stages due to the neutrino emission. 
The mechanism we consider is the subsequent collisions of the produced neutrinos with the outer 
shells of the star. We find that this mechanism can indeed slow down the star rotation but only in 
the first tens of seconds of the core formation, which is when the appropriate conditions of flux and 
collision rate are met. We find that this mechanism can extract less than 1% of the star angular 
momentum, a result which is much less than previously estimated by other authors. 
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I. INTRODUCTION 

Neutrinos play a significant role in the evolution of 
various astronomical objects. They should carry away 
almost all of the gravitational energy lost in the collapse 
of a massive star during a supernova explosion [l[ , phe- 
nomenon that was confirmed by the detection of neu- 
trinos from supernova SN1987A, in the available exper- 
iments at the time, namely Kamiokande-II @, 1MB Q, 
Baksan [|, and LSD @|. Neutrinos are also important 
in the cooling of the neutron star formed at the center 
of a supernova explosion 0, 0]. Neutrinos, because of 
their weak interaction with matter, provide important 
signals of the inner parts of stars like our Sun Q and 
of the cores in stellar collapses Q . Asymmetric neutrino 
emission may also be responsible for the large peculiar 
velocities observed in pulsars [Hi El ■ It has also been 
proposed a long time ago that neutrino emission can slow 
down the spinning of neutron stars [12|,|13|. Here we want 
to revisit the latter idea, trying to refine the estimate by 
using our current knowledge of the physics of neutrinos 
and their emission in neutron stars. 

The majority of neutron stars are known to have large 
angular velocities, and in the case of radio pulsars one 
can directly measure their speed of rotation. It is also ob- 
served that, on average, their rotation tends to slow down 
with time, a phenomenon that is explained by emission 
of electromagnetic waves or, in some conditions, by the 
emission of gravitational waves or other processes [l4| . 
This should be the case during most of the life of the 
neutron star. However, at the early stages, during the 
collapse and formation of the star core (in a time scale 
of 10 s), which is when the flux is intense enough and 
the mean free path is comparable to the star size [l5| . 
the spin-down of this protoneutron star (PNS) can be 
influenced by the collisions of the neutrinos with the star 
matter as they escape. It is interesting to notice that 
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the opposite effect, i.e. the acceleration of a neutron 
star rotation, has also been proposed 16] for the case of 
neutrinos interacting with strong toroidal magnetic fields 
inside stars. 

Note that neutrino emission can also decrease the an- 
gular momentum of the star due to the change in the 
gravitational mass of the star [l7l - fl9j . According to the 
estimate of Ref. [l9| , a star can lose at most 40% of its ini- 
tial angular momentum. Other proposed mechanisms of 
angular momentum loss at early stages of the PNS evolu- 
tion include viscous processes [201 ] . transfer of rotational 
energy into the ener gy o f the supernova explosion 12 1|. 
magnetic PNS winds [22| , and propeller mechanisms [23| . 
However, it was found in Ref. [24[ that none of these the- 
oretical explanations can robustly spin-down a PNS from 
about several ms to the observed periods of rotation of 
young pulsars. We must add that the loss of angular mo- 
mentum by neutrino emission may also occur if there is 
anisotropy at the neutrino production points, which may 
be the case if the star matter has significant polarization 
due to the star magnetic field [25j]. However, we do not 
consider this effect here. In our work we are interested in 
the case in which the production is rather isotropic, but 
neutrinos travel a sizable distance in the star and sub- 
sequently scatter with matter in outer shells where the 
transverse velocity of the medium is larger. 

In this work we revise the previous estimates for the 
spin-down of PNS by neutrino emission [l^, EH , where it 
was stated that this mechanism can even possibly stop 
the star rotation. This result denotes a very dramatic 
effect, which we want to study using a more detailed 
work, but still within analytical models in order to study 
the sensitivity of our results to different parameters of a 
generic star. Our estimates show that the effect is much 
weaker than previous estimates - less than 1% reduction 
of the angular speed. We also check that most refine- 
ments in our treatment point to further reduction, not 
enhancement, of the effect. 

In Sec. [H] we present our conceptual treatment and cal- 
culation of the spin-down of PNS due to neutrino emis- 
sion. In Sec. Mil we summarize our results and state our 
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conclusions. 



II. MODEL AND CALCULATIONS 

In this section we formulate the phenomenon of spin- 
down of a forming neutron star due to neutrino emission 
during the core collapse. The spin-down mechanism con- 
sists in the fact that neutrinos which escape from the star 
are produced in regions around the neutrinosphere, where 
the transverse velocity of matter is relatively lower, and 
as they propagate to outer regions, some of them will col- 
lide with matter moving with larger transverse velocities, 
thus absorbing transverse momentum and causing a spin 
down of the medium. In this sense, the trajectory of an 
outgoing neutrino should bend as it propagates through 
the rotating medium, due to collisions. For this effect 
to be of any significance, the mean free path should be 
less (but not much less) than the star radius. Otherwise, 
if the mean free path is much larger, there will be too 
few collisions, while if it is much shorter, the difference 
in transverse velocities from the emission to the collision 
points will be too small. 

As cited above, spin-down due to neutrino emission 
was already proposed in the past. However, with the 
knowledge available today we can include more details 
in the treatment, namely (i) take proper account of the 
weak interactions in the collisions, (ii) use a phenomeno- 
logical matter composition and density profile in the star 
instead of a uniform medium, (iii) take into account the 
opacity of inner parts of the star where neutrinos ther- 
malize, so for the spin-down effect neutrinos are emitted 
only from a neutrinosphere instead of from the center of 
the star, and (iv) use thermal spectra for each neutrino 
species instead of a single monochromatic emission. 

We organize the calculation starting from the neutrino 
production at the neutrinosphere, followed by the indi- 
vidual neutrino collisions with the rotating medium fur- 
ther outside, and finally, by adding these collisions all 
around the star, to get the result. At each step we state 
the model of the situation and the approximations used. 

Neutrinos are produced everywhere around the star, 
but at high temperatures (energies) they suffer much 
scattering and absorption. Therefore, those which man- 
age to escape are not produced at the star center but at 
the neutrinosphere, or surface of last scattering (26|. The 
definition of the neutrinosphere is statistical and depends 
on the medium density as well as the neutrino species and 
energy. Here we will simply define it as a spherical shell 
of radius R ns , different for every neutrino species and 
energy, given by one mean free path less than the star 
radius, 
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where R is the star radius and a v is the (energy- 
dependent) total cross section for neutrino scattering in 
the star medium with nucleon number density n. Since 



n depends on the position, Eq. (|2.1[) is really an equation 
for R ns that we solve in each case. We consider a PNS 
medium where a fraction Y n ~ 0.9 of the nucleons are 
neutrons and a fraction Y p ~ 0.1 are protons [(|. 

We neglect scattering with electrons, because they are 
in equal number to protons but their cross section is an 
order of magnitude smaller. We also neglect Pauli block- 
ing or nucleon correlations. These effects should be more 
important at later stages of the star evolution, when neu- 
trino energies are lower; however, the spin-down caused 
by neutrinos is negligible then. In any case, correlations 
tend to reduce the cross sections [27|, pointing further 
into the direction that previous calculations of the PNS 
spin-down by neutrinos were overestimated. 

The total cross sections (and thus neutrinospheres) dif- 
fer for three species of neutrinos: electron neutrinos v e , 
electron antineutrinos v ei and all other v x . The neutrino- 
nucleon cross sections, for neutrinos with energies 1 — 100 
MeV, are approximately (see pp. 160-167 in Ref. [l[) 

C {v e n)inclastic —a(p e -p) inelastic 

o-(^)dastic =2.6 x 10~ 42 | Ev ) cm 2 , 
v Elastic ^10 MeV/ 

(all v species), 

crMeiastic =2.1 x lO" 42 ( — ] cm 2 , 
v ^da^ic ^10 MeV/ 



(all v species) 



(2.2) 



The total cross sections in the star medium for the species 
in question are then 

<?u e = Y n [a(v e n)i ne i, + a{vn) e \.) + Y p a(vp) Q \., 
c Pc = Y p [a(i>eP)inel. + v(vp)eh] + Y n <j{vn) c \., 
v Vx = Y n a(un) e i. + Y p a(vp) c \., (2.3) 

which result in the hierarchy er„ c > er Pc > a Vi . Conse- 
quently, the sizes of the respective neutrinospheres follow 
the same order and, since the star temperature is higher 
further inside, the average energies also follow a hierar- 
chy [2g| : 



(E Ve ) ~ 10 MeV, 



20 MeV. 



15 MeV, 



(2.4) 



Concerning the neutrino spectra, we take two alternative 
approaches. In our first, simpler approach, we consider 
purely monoenergetic neutrinos for each species, as given 
in Eq. (JTTJ) . Each species Vi has then a definite neut- 
rinosphere. In our second, more refined approach, we 
consider thermal (Fcrmi-Dirac) energy distributions for 
each neutrino species [29j], with temperatures according 
to the energy averages given in Eq. pi]). In this case, the 
neutrinospheres are continuously distributed. 

We also checked for the possibility that neutrinos could 
experience flavor oscillations in matter while propagating 
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from the neutrinosphere towards the star surface, but 
found this effect to be irrelevant for the spin-down. Neu- 
trino resonant conversion at energies near 10 MeV is im- 
portant in the expanding envelope [30| where the matter 
density is (10 26 — 10 27 ) cm~ 3 , but in our region of interest 
of the star densities are about 10 35 cm -3 (see Ref. 
where flavor oscillations are suppressed. Oscillations can 
happen, but in a rather thin layer close to the neutron 
star surface, and so their effect on the spin-down is neg- 
ligible. 

Another group of flavor changing processes, which can 
influence the angular momentum transfer by neutrinos, 
is the neutrino flavor conversion due to the v — v scat- 
tering [3ll Hij. This effect was shown to significantly 
change the initial flavor content of neutrinos in a dense 
neutrino flux, corresponding to the neutrino luminosity 
C u > 10 51 erg/s. It was however found in Ref. [12] that 
the significant transition probability due to v — v colli- 
sions is archived at the distances (70 — 80) km from the 
neutrinosphere surface, i.e. it happens in the envelope 
of a star. On the contrary, we study the spin-down of 
PNS due to the collisions with background matter in the 
core of PNS outside the neutrinosphere at the distances 
< 20km from the star center (see Fig. [2] below). 

Note that background matter can also influence the 
collective neutrino flavor transformations. In order to 
have some resonance effects in these flavor changing 
processes, the neutrino density n„(r) = C l ,/4:TTr 2 (E u }, 
where (E v ) ~ 10 MeV is the typical neutrino energy and 
C v ~ 10 52 erg/s, should be comparable with the elec- 
trons number density n e . In our calculations we suggest 
that n e « 0.1n„ ~ 10 34 cm -3 inside the core of PNS at 
r < (15 — 20) km. The number density of neutrinos at 
the same distance is n v ~ 10 32 cm -3 , which is 2 orders 
of magnitude less than the electron density. It means 
that the background matter is unlikely to generate any 
resonance effects in our case. 

Now we should address the density profile of the star. 
The actual density profile strongly depends on the equa- 
tion of state of the nuclear matter, so these density dis- 
tributions are generally not well known [33j. Results of 
numerical simulations for the density profiles [34| can be 
approximated by the following expression: 

n(r) =n c (l- cxp > ( 2 -5) 

where n c is the central density and e is a phenomenolog- 
ical parameter. The e = case corresponds to the well 
known Tolman VII model [35j |. An analogous expression 
for the density profile was used in the study of nonradial 
oscillations of a neutron star [36j]. In Fig. QJa) we show 
the behavior of the density for various values of the pa- 
rameter e. In order not to obscure the sensitivity of our 
results to different parts of the calculation, yet with the 
risk of not being realistic in specific cases, we will use 
these analytical density profiles instead of full numerical 
profiles. As shown in Fig. [Ha), negative e implies a more 
flat profile with a fast drop at the surface, while positive 



values imply a gradual decrease of the density from the 
center to the surface. 

From Eqs. (f2~Tj) . ([2~B1 and (l2"3j) we can get the radii 
of the neutrinospheres as functions of the parameter e, 
for a given neutrino energy. In our analysis we use the 
values n c = (4.7 — 6.6) x 10 38 cm~ 3 and R = 15 km as 
well as n c = (2.0 - 2.7) x 10 38 cm" 3 and R = 20 km for 
the central density and the star radius. The total mass 
of the star calculated on the basis of Eq. (|2.5p has the 
form, 

F M (e)= dxx 2 (l~x 2 )e- ex2 . (2.6) 
Jo 

For the chosen densities and radii Eq. (|2.6[) gives one the 
mass of PNS in the range (1.4 — 2.0)M Q , depending on 
the parameter e. Although the mass of a neutron star 
strongly depends on the equation of state of the neu- 
tron star matter, it is unlikely bigger than 2M Q [37|. In 
Fig. [T^b) we present the e dependence of the function 
Fm{^)i which is proportional to the PNS mass. 

Now, let us address the issue of emission directions. 
From the neutrinosphere, neutrinos are emitted outwards 
in all directions. However, as a simplified model for the 
spin-down effect, we approximate the emission as purely 
radial, where "radial" is meant in the frame comoving 
with the local matter at the neutrino production point. 
In reality, of course, neutrinos are produced in all di- 
rections and rigorously one should consider the whole 
hemisphere of outgoing directions at every point. How- 
ever, those particles emitted more towards the tangential 
velocity will transfer less momentum to the star under 
subsequent scattering, while those emitted against the 
tangential velocity should transfer more momentum. In 
order to test the validity of the radial emission model, 
we performed the calculations in the case of neutrinos 
emitted in all directions within the equatorial plane and 
checked that indeed collisions with purely radially emit- 
ted neutrinos represent the full spin-down effect up to a 
geometrical factor near unity. Since there is no need for 
precision higher than a few tens of a percent, in what 
follows we will present the estimates within the model of 
radially produced neutrinos. 

The process is then calculated by considering that neu- 
trinos are emitted from their corresponding neutrino- 
sphere and subsequently collide with the star medium 
on their way out, thus taking away part of the angular 
momentum. 

Now let us describe the subsequent elastic collision 
of the outgoing neutrino with the star matter, which is 
where the spin down effect takes place. The effective 
weak interaction for elastic neutrino scattering with a 
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(a) (b) 



FIG. 1. (a) The radial dependence of the density for various values of the parameter e. The solid line corresponds to the 
Tolman VII model (e = 0), the dashed line is built for e = 1, and the dash-dotted for e = — 1. (b) The function Fm(s), given 
in Eq. 1|2.6[> . which defines the star mass dependence on the density profile parameter e. 



TABLE I. The values of the coefficients g L ,R in Eq. (12. 7|) for 
various channels of neutrino elastic scattering. v x stands for 
x 7^ e while v stands for all lepton flavors. 



No. 


Reactions 


9L 


9R 


I 


v £ e - 


-> v<,e 


0.73 


0.23 


2 


v m e - 


» v x e 


-0.27 


0.23 


3 


v £ e - 


* v e e 


0.23 


0.73 


4 


v x e - 


-> v x e 


0.23 


-0.27 


5 


vp - 


■¥ vp 


0.27 


-0.23 


6 


vn - 


■> vn 


-0.50 


0.00 


7 


vp - 


->• vp 


-0.23 


0.27 


8 


vn - 


-> vn 


0.00 


-0.50 



background fermion / has the form [36 



G 



V2 J 

x P(fc 2 ) 7 "(l - 7 £ >(fc 1 ), 



(2.7) 



where v(ki^) are the initial and final neutrino spinors 
with momenta fci j2 = (wi, 2)^1,2), and /(pi.2) are the 
respective spinors of the fermions in the medium, with 
momenta px >2 — (-E1.2, Pi, 2) respectively. The coefficients 
9l,r depend on the neutrino scattering channel and are 
listed in Table Q] 

From the effective interaction of Eq. (|2.7p , we can build 
the cross section, 



da(vf -> vf) =-^r^5 A (p2 + k 2 -p\- k\) 



64tt 



IMI 



-d 3 k 2 d 3 p 2 , 



(2. 



(fci - pi) uj 2 E 2 

where the matrix element squared for unpolarized scat- 



tering derived from Eq. (|2.7|) is 

\M\ 2 =128G|[ 5 |(p 1 • h) 2 + g 2 R ( Pl ■ k 2 f 

-9L9Rm 2 f {k 1 -k 2 )} ) (2.9) 

with mf being the mass of the fermion /. For radially 
emitted neutrinos, and taking this radial direction at each 
collision point as our reference, the incoming fermion 
from the medium will be perpendicular and the outgoing 
neutrino will have an arbitrary direction parametrized by 
the relative polar and azimuthal angles 62 and 4> 2 , such 
that 

ki =wi(0,0,l), pi =£117(0,1,0), 

k 2 (sin 02 cos (fo, sin 02 sin 02) cos #2)- (2.10) 

Here Vf is the tangential velocity of the star at the in- 
teraction point. We are assuming that the thermal ve- 
locities of the fermions in the medium average out con- 
cerning this effect, and the only effective velocity is the 
average drift of the medium due to rotation. With these 
conventions, the dot products of interest are 

(pi • ki) =£Wi , (ki ■ k 2 ) = UiL) 2 (l - cos0 2 ), 

(pi • k 2 ) =E x u) 2 {l - v f sin 2 sin</> 2 )- ( 2 -H) 

We can integrate Eq. (I2.8[) over P2 to get the differential 
cross section, 



da _ \M\ 2 wjf 



dn 2 64tt 2 u\E 2 ' 



(2.12) 



where the energy of the outgoing neutrino is given, by 
energy-momentum conservation, in terms of the scatter- 
ing angles 



W 2 = 



E\{Y — Vf sin0 2 sin 2 ) -I- wi(l — cos 02) 



(2.13) 
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Now we can calculate the rate of transverse momentum 
transferred to the star at that collision point. Let J(r) 
be the (radial) neutrino flux incoming to a given collision 
point; then J(r)(da / 'df^) df^ is the rate of outgoing neu- 
trinos within . Since the component of momentum 
of each of these outgoing neutrinos is u>2 sin 82 sin cf>2 , then 
the total momentum per unit time transferred to the 
star due to the outgoing neutrinos ejected from a given 
collision point is 



(kj,) = J uj 2 sin #2 sin 4> 2 J(r) 



da 



da 



'2, 



(2.14) 



where dSl2 = sin^d^d^- The flux J(r) at the collision 
point is related to the neutrino flux at the surface of the 
star, J , by the relation J(r) = J x (R/r) 2 , where R is 
the star radius. Using Eq. (|2.14p we can now calculate 
the total rate of angular momentum transferred to the 
star as a whole by just summing over all collision points 
that lie outside the corresponding neutrinosphere, 



L z = I (fc )rsintf n f (r) d 3 

'r>R ns 



(2.15) 



where n / (r) is the local number density of target fermions 
in the medium, and $ is the polar angle (colatitude) of 
the star at the collision point. The full result for L z 
is the sum of this type of calculation, repeated for each 
neutrino species. 

The rate of angular momentum transfer to the star as 
given in Eq. (I2.15[) requires the computation of a four-fold 
integral which has to be evaluated numerically. Never- 
theless, within the following approximations one can cal- 
culate an analytical expression for L z . First, assuming a 
star radius R — 20 km and angular velocity = 10 3 s" 1 , 
its equatorial linear velocity (in units of c) is ~ 0.1. Thus 
we can treat Vf in Eq. (|2.13[) as a small parameter. In 
addition, typical energies for the emitted neutrinos are 
lo\ ~ 10 MeV, while the incoming fermion in the medium, 
which is nonrelativistic, has an energy E\ near the nu- 
cleon mass. Therefore the ratio u>\/E\ is also small. Con- 
sequently, we use the approximation, 



U> 2 «Wl(l + £)) 



(2.16) 



where 



£=Vf sin 9 2 sin ^2 — ^(1 



- COS 62) < 1. 



For the matter velocity at the collision point we use 

^* R'ns 



Vf — vq sin •& 



R 



(2.17) 



where vo is the equatorial velocity of the neutron star. 
Here Vf corresponds to the relative transverse velocity of 
the medium at the collision point (radius r) with respect 
to the velocity of the medium at the neutrino emission 
point (radius R ns )- 



Using Eq. (|2 .9[) and Table Q] we can determine the 
square of the matrix element, for example for the spe- 
cific reaction v e n — > v e n, which is \M\ 2 — 32GpE 2 Lu 2 . 
Consequently, using Eqs. (|2.15[) and (|2.16|) one gets the 
following approximation for the rate L z due to this type 
of collision: 



L z ^2GlujfRJ v / r(r - R ns )dr 

x f n„(r,^)sin 3 ?9di?, 
Jo 



(2.18) 



where n n (r, 1?) is the neutron density at the correspond- 
ing collision point. 

It is convenient to present the final result as the ratio 
between the rate of angular momentum loss and the ini- 
tial angular momentum of the star Lq = where its 
moment of inertia is given by 



(2.19) 



X= — m n [ drr 4 n(r), 
3 Jo 



and where m n is the neutron mass. Based on Eqs. (|2.5 
and (|2.18l) we find for this ratio: 



L z _ G\L V E 2 Fi(e) 
T ~4ir 2 R 2 m n F (e) 



1.0s" 1 x 



E„ 



R 



10 MeV / V 1052 er g/ S 



lOkmy F (e)' 



(2.20) 



where C v is the neutrino luminosity and where we have 
defined the following integrals: 



F (e) = / dxx 4 (l-x 2 )e- ex , 
Jo 

Fi(e)=/ dxx(x~x ){l-x 2 )e- tx2 

J XO 



(2.21) 



with xq = R ns /R. 

We have just considered the reaction v e n — > v e n. The 
contributions of all other neutrino species, including col- 
lisions with protons, are treated in a similar way and 
added to the result. The average neutrino luminosity just 
after the neutronization stage [28[ is C v ~ 10 52 erg/s. 
This large neutrino luminosity lasts for a few seconds, 
mainly during the Kelvin-Helmholtz stage, with all neu- 
trino types, v e41T and v e .^.T, having almost equal lu- 
minosities, and therefore giving similar contributions to 
Eq. (I2~20l) . 

A refinement of this neutrino emission model considers 
the three neutrino species not monoenergetic but with 
thermal distributions, each with distinct temperatures 
and chemical potentials. Therefore we can choose the 
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Fermi-Dirac distribution in the form [2i 



dN 
dE~ 



El 



F(t) v )T* exp{E v /T v - rh,) 
dx- 



V 



e x-r, v + I 



(2.22) 



where we have the relation (E v ) /T v « 3. 1514+0. 125077^+ 
0.042977 2 , + • • • between the temperature of the neutrino 
gas T v and the mean neutrino energy (E v ) defined in 
Eq. ((11]). The typical values of the chemical potentials 
are (see Ref. |29j |) r) v<s ~ 2, rj 9c ss 3, and r]„ x s» 1. 

On the basis of Eqs. ([2~2Tty - (j2~2"2"]) we get the averaged 
ratio as 




«1.0s- 



10 MeV / 



R 



F( v , 



10 52 erg/s J Vl0km / 
1 f°° y 5 dy F^cy) 
n e»-^+l F (e) 



(2.23) 



where we are including the energy dependence of the 
function Fi(e) in Eq. (|2.21[) in the form of the dimen- 
sionless parameter y — (J5 w /10MeV), since the size of 
the neutrinosphere depends on the neutrino energy [see 
Eqs. and 1|23J)]. 

In Fig. [5] we present the total angular momentum car- 
ried away by all neutrino species, L z — ^2 i L z (vi)At, 
where At = 10 s , in units of the initial angular momen- 
tum Lq, versus the mass of PNS. Here At is the time 
interval during which the luminosity stays at the high 
value L v ~ 10 52 erg/s. Note that the PNS mass varia- 
tion on this plot results from its epsilon dependence [see 
Eq. (EU) and Fig.[TJb)]. 

As one can see from the figures, a larger fraction of 
the angular momentum can be carried away by neutri- 
nos for a star with bigger radius. We can also notice 
that with the enhancement of the central density the ef- 
fect increases. On one hand, the neutrinosphere radius 
grows with the enhancement of the central density (|2.1I) . 
Thus neutrinos have less opportunity to collide with ro- 
tating matter and the relative velocities from production 
to collision are also smaller. On the other hand, PNS 
with equal masses and different central densities corre- 
spond to a different e parameter in Eq. (|2.6|) . For ex- 
ample, n c = 4.7 x 10 38 cm -3 corresponds to e = 1 and 
n c = 6.6 x 10 38 cm" 3 to e = 2 for M = 1.4M Q . The big- 
ger value of e implies a more steep density profile near the 
central region and, hence, a bigger concentration of PNS 
mass there. It means that effectively a neutrinosphere 
should decrease at higher e. The results of numerical 
simulations presented in Fig. [5] shows that the latter ef- 
fect is more significant. 

One can also see in Fig. [5] that the difference between 
the models of monoenergetic and thermally distributed 



neutrinos is less than several percent. Thus, a monoen- 
ergetic neutrino model is still a good approximation to 
describe the spin-down. 

Including general relativity corrections to the moment 
of inertia of the neutron star, the values of L z /Lq become 
slightly smaller than those we obtained. For example, if 
we consider a neutron star with M = 2M Q , the reduc- 
tion factors within the Tolman VII model (e = 0) will be 
« 0.76 for R = 15km and 0.82 for R = 20km [Hj]. 
As a result, we find that, by the mechanism we have 
considered, and within the Tolman VII model as bench- 
mark (e = 0) for a R = 15 km star, neutrinos could carry 
away up to 2.4 x 10 of the initial angular momentum of 
the PNS [compare this value to the corresponding result 
shown in Fig. [TJa) , which does not include the general 
relavivity corrections] . 

For other density profiles, i.e., corresponding to M = 
1.4M Q (e = 2), neutrinos can carry away 1.7 x 10~ 3 of 
the initial angular momentum for a star with R = 20 km 
and n c — 2.8 x 10 38 cm~ 3 . This result is significantly less 
than predicted in the previous works [T3 . Il3j ] . 

What follows at later stages of the neutron star evolu- 
tion is that the neutrinospheres shrink as the star cools 
down, because the mean free path grows beyond the star 
radius. This is mainly due to the Pauli blocking reduc- 
tion of the neutrino scattering cross section. At the same 
time, the neutrino flux drops down dramatically from its 
initial values. The expression in Eq. (12.201) is then valid 
only at the initial stages, which is the only time when 
the spin-down due to neutrino emission can be of any 
significance. 



III. SUMMARY AND CONCLUSIONS 

We have studied the spin-down of a forming neutron 
star due to neutrino emission. The spin-down by neu- 
trino emission is significant only for the first few seconds 
of the neutron star evolution, when the size of the neut- 
rinosphere is less than (but comparable to) the radius of 
the star and the neutrino flux is still large. These condi- 
tions allow for neutrinos to have a large enough collision 
rate. At later stages the star becomes almost transpar- 
ent for neutrinos and the neutrino flux is too small for 
this effect to be of any significance, so other spin-down 
mechanisms take over. 

We model the phenomenon by considering neutrinos of 
type v e , i> e , and v x (where x stands for all other), which 
have different interaction with the medium, each pro- 
duced at their corresponding neutrinosphere and moving 
radially outwards, subsequently colliding with the star 
matter in the outer shells, where the transverse velocity 
of the medium is larger than at the production point, 
thus causing a slowing down of the star. 

An arguable part of the model is the consideration 
of purely radially moving neutrinos at the production 
points, but we have checked analytically in some simpli- 
fied cases that the emission in all directions causes an 
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FIG. 2. The fraction of the total angular momentum carried away by all neutrinos types versus the PNS mass, which depends on 
the parameter e according to Eq. (|2.6|l . for various central densities. The dashed lines correspond to the simple model (|2.20|) with 
monoenergetic neutrinos with energies given in Eq. ([II]). The solid lines correspond to the refined model (|2.23jl with thermally 
distributed neutrinos with average energies (temperatures) given in Eq. ((II]). (a) The fraction for a star with R = 15 km. The 
upper curves correspond to n c = 6.6 x 10 38 cm -3 and the lower ones to n c — 4.7 x 10 38 cm -3 , (b) The fraction for a star with 
R — 20km. The upper curves correspond to n c = 2.7 x 10 38 cm -3 and the lower ones to n c = 2.0 x 10 38 cm~ 3 . 



average effect not much different than the purely radial 
case. 

Another simplification is that Pauli blocking and other 
nucleon correlations are not taken into account, presum- 
ing that they are more important at later stages, when 
neutrino flux and energies are lower and this spin-down 
mechanism is negligible. 

Another important part of the model concerns the den- 
sity profile, where we use a phenomenological analytical 
expression instead of more realistic numerical profiles, in 
order to study the sensitivity to it and have a better 
comparison to previous estimates to the spin-down phe- 
nomenon. We find that the density distribution, n(r), 
significantly affects the results for the rate of angular 
momentum loss (see Fig. [I}. Very few analytical for- 
mulas for n(r) models are known, and the majority of 
the density profiles are available from numerical simula- 
tions. We have chosen an exponential density profile (|2.5[) 
which depends on a phenomenological parameter e to fit 
the results of the numerical simulations [sij . For e = 
this formulation reproduces the well known Tolman VII 
model. Although this model cannot capture the fine de- 
tails of individual cases, it reproduces the main features 
of realistic stars. However, it is unable to explain without 
fine tuning the enormous range of matter densities, from 
~ 10 6 g/cm 3 in the crust to > 10 14 g/cm 3 in the center. 
Therefore we have introduced an additional exponential 
factor in the model for n(r). 

As shown in Fig. HJa), for positive e the density goes 
gradually down, while for negative e it stays high almost 
all over the star, falling sharply at the surface. Consis- 
tently, Fig. [2] shows that the spin-down effect is bigger 
for PNS with smaller mass at the fixed central density. 
Taking into account Fig. QJb) we obtain that spin-down 
should increase with the enhancement of e because a 



higher and flatter density implies a larger neutrinosphere 
radius, thus decreasing the relative velocity of the star at 
the collision point with respect to the point of produc- 
tion. On the contrary, if neutron stars would feature a 
more gradual density descent, the average distance from 
the neutrinospheres to the collision points will increase, 
causing a stronger spin-down effect. 

We have found that, for some density profiles corre- 
sponding to e = 2 in the case of a PNS with R = 20 km 
and central density n c = 2.8 x 10 38 cm~ 3 , which gives 
one M — 1.4Mq, neutrinos can carry away up to about 
1.7x 10~ 3 of the initial angular momentum of the neutron 
star, provided that an average neutrino luminosity near 
C v = 10 52 erg/s lasts for about 10 s during the Kelvin- 
Helmholtz stage of the neutron star evolution. The re- 
sults of our calculations for other central densities and 
radii are presented in Fig. [2] 

We should recall that previous estimates [HI, [l3[ gave 
larger values of L z , and it was even predicted that the 
rotation of a neutron star could be stopped by neutrino 
emission (or equivalently, a reduction of the angular ve- 
locity by more than an order of magnitude) . After using 
a variable density profile, our better knowledge of elec- 
troweak interactions today, and the most recent data on 
the neutrino flux, we find that the results are clearly not 
as dramatic. 

Let us examine the importance of various factors which 
were not accounted for in the previous calculations of the 
PNS spin down [l2l . [l3j . For example, we can discuss 
the situation when all neutrinos are monoenergetic, with 
E v ~ 10 MeV and C u ~ 10 52 erg/s, and are emitted from 
the PNS center, whereas the density profile is given by 
Eq. (|2.5p . The considered case is equivalent to the neut- 
rinosphere with zero radius: xq = in Eq. (12.211) . Using 
Eq. (|2~2H|) for PNS with R ~ (15-20) km, we obtain that 
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L z At ~ Lq, where At ~ several seconds. Therefore we 
obtain that the rotation of PNS can be significantly re- 
duced by the neutrino emission provided all the particles 
are emitted in the center of PNS, which reproduces the 
results of Refs. [IH, It means that the concept of the 
neutrinosphere is the most important in our calculations. 

It can be also shown that various density distribution 
profiles can change the PNS spin-down, however not so 
dramatically. Concerning other corrections, the inclusion 
of Pauli blocking or nucleon correlations would tend to 
make the effect even smaller; different proportions of neu- 
trino species or neutrino oscillations only induce minor 
or negligible changes; using monoenergetic or thermally 
distributed neutrinos also result in minor differences in 
the spin-down effect. 

Concerning the experimental observation of this effect, 



unfortunately there is quite limited information about 
the initial angular velocities of neutron stars. An ef- 
fort to infer initial angular velocities of PNS has been 
made [3{|, where it was revealed that there should be a 
large uncertainty in the results. 
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